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The G1/S transition is a crucial decision point in the cell cycle. At G1/S, there is an abrupt 
switch from a state of high cyclin-dependent kinases (CDK) inhibitor (CKI) levels and 
low S-phase CDK activity to a state of high S-phase CDK activity and degraded CKI. 
In budding yeast, this transition is triggered by phosphorylation of the Cdk1 inhibitor 
Sid at multiple sites by G1-phase CDK (Clnl ,2-Cdk1 ) and S-phase CDK (Clb5,6-Cdk1 ) 
complexes. Using mathematical modeling we demonstrate that the mechanistic basis 
for the abruptness of the G1/S transition is the highly specific phosphorylation of Sicl 
by S-phase CDK complex. This switch is generated by a double-negative feedback loop in 
which S-CDK1 phosphorylates Sicl , thus targeting it for destruction, and thereby liberating 
further S-CDK1 from the inhibitory Sic1-S-CDK1 complex. Our model predicts that the 
abruptness of the switch depends upon a strong binding affinity within the Sic1-S-CDK 
inhibitory complex. In vitro phosphorylation analysis using purified yeast proteins revealed 
that free Clb5-Cdk1 can create positive feedback by phosphorylating Sicl that is bound in 
the inhibitory complex, and that Sid inhibits Clb5-Cdk1 with a sub-nanomolar inhibition 
constant. Our model also predicts that if the G1-phase CDK complex is too efficient 
at targeting Sid for destruction, then G1/S becomes a smooth and readily reversible 
transition. We propose that the optimal role for the G1-phase CDK in the switch would 
not be to act as a kinase activity directly responsible for abrupt degradation of CKI, but 
rather to act as a priming signal that initiates a positive feedback loop driven by emerging 
free S-phase CDK. 
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INTRODUCTION 

The major cell cycle transitions are triggered by phosphorylation 
or dephosphorylation of the targets of cyclin-dependent kinases 
(CDKs) (Morgan, 2007). At the transition from Gl- to S-phase, 
the S-phase promoting signal is generated by downregulation 
of CDK inhibitors (the CKIs) that suppress the activity of the 
S-phase CDK (Peter, 1997). In budding yeast, the CKI acting 
at the Gl/S transition, the protein Sicl, is simultaneously both 
an inhibitor and a substrate of Cdkl. At the onset of S-phase 
Sicl is phosphorylated by Cdkl, thereby generating two "diphos- 
phodegrons" containing two properly spaced phosphates that are 
recognized by the Cdc4-SCF ubiquitin ligase (Hao et al., 2007). 
Ubiquitination by Cdc4-SCF targets Sicl for proteolysis via the 
proteasome pathway (Verma et al, 1997b). 

Sicl was first identified as a substrate of Cdkl by immunopre- 
cipitation (Reed et al, 1985; Wittenberg and Reed, 1988). Later, it 
was shown to be an inhibitor of Clb-Cdkl complexes that reg- 
ulates cell cycle progression at the M/Gl and Gl/S transitions 
(Mendenhall, 1993; Donovan et al., 1994). Transcription of SIC1 
starts in late mitosis and its protein levels increase until the end 
of Gl phase, followed by a rapid turnover at the Gl/S transi- 
tion (Schwob et al., 1994; Verma et al., 1997b). The signal for the 
degradation was proposed to be the multisite phosphorylation 
of Sicl by Gl-specific Clnl,2-Cdkl complexes. This is possi- 
ble because, in contrast to S- and M-phase specific Clb-Cdkls, 



Cln-Cdkls are not inhibited by Sicl (Verma et al., 1997a,b; 
Nash et al, 2001). The destruction of Sicl was shown to release 
the Clb5,6-Cdkl complexes required to initiate DNA replication 
(Schwob and Nasmyth, 1993; Schwob et al., 1994; Schneider et al., 
1996). 

Recently, we showed that the Cdkl -dependent phosphoryla- 
tion of Sicl at multiple sites occurs in a specific sequence that 
forms a semi-processive phosphorylation cascade. This cascade 
is mediated by Cksl, the phosphate binding adaptor subunit of 
Cdkl (Koivomagi et al., 2011a). The S-phase promoting com- 
plex Clb5-Cdkl was shown to be able to efficiently phosphorylate 
all the crucial sites on Sicl that were required for generation of 
phosphodegrons. On the other hand, the Gl-specific Cln2-Cdkl 
complex, while still capable of phosphorylating Sicl, was much 
less efficient at transferring phosphates to the critical phosphode- 
gron residues at the end of the cascade. Importantly, Cln-Cdkl 
complexes alone were not able to cause the degradation of Sicl 
in vivo. This led us to propose that positive feedback, arising from 
the emerging free Clb5,6-Cdkl complexes, is responsible for the 
degradation of Sicl, while the Clnl,2-Cdkl complexes play an 
auxiliary role as a primer of the phosphorylation cascade. 

The relative effectiveness of either Clns or Clbs as the driv- 
ing force of the Gl/S switch depends on several factors, including 
relative rates of Sicl phosphorylation and the strength of inhi- 
bition of Clb complexes by Sicl. Therefore, to understand the 
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switch, it is crucial to combine quantitative biochemical analy- 
sis of these parameters with a modeling approach to evaluate the 
optimal parameters with respect to the dynamic qualities of the 
switch. To function as a sharply defined transition point between 
different cellular states, the CKI- and CDK-controlled Gl/S tran- 
sition should possess both ultrasensitivity and high degree of 
irreversibility. An ultrasensitive signal-response filters noise by 
setting a sharp threshold for switch activation and being non- 
responsive to small signal fluctuations. An abrupt transition is 
also important for coherence in the activation of S-phase pro- 
cesses and regulatory loops (Skotheim et al, 2008). For example, 
a gradual transition could create a window in which enough Cdkl 
activity exists to fire some origins of replication but not enough 
activity to fully inactivate all the origins and thereby prevent 
genotoxic reinitiation. A high degree of irreversibility would also 
ensure that transient decreases of stimulus would not reverse the 
switch and interrupt processes that have already initiated. 

A sharp transition threshold and irreversibility can be achieved 
by creating bistability in the switching system (Novak et al., 
2007; Domingo-Sananes and Novak, 2010; Zhang et al., 2011). 
Bistable switches have been described for cell cycle transitions 
(Ferrell et al, 2009; Trunnell et al., 201 1), and the double-negative 
feedback system of the CKI/CDK module governing the Gl/S 
transition is a potential source of bistability. Alternatively, delayed 
accumulation of multiply phosphorylated forms of Sicl has been 
proposed as a distinct mechanism to create ultrasensitivity at the 
Gl/S transition (Nash et al, 2001), and also for other switches 
(Kim and Ferrell, 2007; Takahashi and Pryciak, 2008; Harvey 
et al, 2011; Trunnell et al, 2011; Lu et al, 2012). Again, it is 
not clear which of the alternative mechanisms provides the opti- 
mal solution in the CKI/CDK system involving two different 
cyclin-CDK complexes and a CKI. 

In the present study we analyzed the optimal range and com- 
binations of some of the kinetic parameters governing the Gl/S 
switch. Although a number of parameters may have an impact 
on the switch, including protein synthesis rates, ubiquitinyla- 
tion rates, proteolysis rates, and transcription factor (SBF/MBF) 
activation rates, in this study we focus on the analysis of param- 
eters that control the mutual counteraction of CDK and CKI. 
We used a minimal model containing both Gl- and S-phase 
CDK complexes together with the CKI, which is an inhibitor 
of the S-phase-CDK. For the model simulations we used our 
experimentally validated range of inhibition and phosphorylation 
parameters and realistic concentrations based on experimental 
data obtained from the literature. The results provide several 
important predictions about the parameter ranges that allow for 
bistability of the system, as well as the realistic dynamic paths 
of the Gl/S switch. The conclusions based on these simulations 
lay out a general framework for understanding the CDK/CKI 
switches in eukaryotes. 

RESULTS 

Sicl INHIBITS Clb5-Cdk1 ST0ICHI0METRICALLY, WITH 
SUBNAN0M0LAR K, 

One of the important parameters controlling the Gl/S switch is 
the strength of the inhibition of Clb5/6-Cdkl complexes by Sicl. 
To quantitatively measure this inhibition we performed an in vitro 



phosphorylation and inhibition assay with purified Sicl and 
cyclin-Cdkl complexes. Figures 1A-C presents the dose-response 
curves showing the effect of different concentrations of Sicl on 
the initial velocity of histone HI phosphorylation by Cln2-Cdkl 
and Clb5-Cdkl. The IC50 value for Clb5-Cdkl was 1.7 nM, which 
places Sicl into the category of tight stoichiometric inhibitors. 
As expected, the Gl complex Cln2-Cdkl was not inhibited by 
Sicl, and an effect was apparent only in the micromolar range 
that corresponds to the Km value for the phosphorylation of a 
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FIGURE 1 | Inhibition of Clb5-Cdk1 and Cln2-Cdk1 by Sicl. (A) 

Autoradiograph of a Clb5-Cdk1 and Cln2-Cdk1 inhibition assay. Sicl is 
extensively phosphorylated by Cln2-Cdk1 whereas no detectable 
phosphorylation is seen with Clb5-Cdk1. (B) Inhibition curves for 
C!b5-Cdk1 and Cln2-Cdk1 are presented as the initial rates of histone H1 
phosphorylation in the presence of different concentrations of Sid. (C) 
Calculated IC50 values from the experiment in panel (A). (D) The values of 
IC50 change in parallel with Clb5-Cdk1 concentrations when both Sicl and 
the kinase are varied in the picomolar or low nanomolar ranges. 
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non-inhibitory version of Sicl (data not shown). This suggests 
that Sicl only inhibits Cln2-Cdkl as a competitive substrate. 

The experiment presented in Figures 1A-C contained 500 u,M 
ATP and 2.5 [iM Histone HI. The K M ATP = 70 [iM, and the value 
for Kjvi.hi is above 10 |xM (our unpublished data). Therefore, the 
obtained value of IC50 is a valid estimate of the apparent inhi- 
bition constant for substrate phosphorylation at saturating and 
nearly physiological ATP concentration. However, the experiment 
presented in Figures 1A-C was performed at low Clb5-Cdkl lev- 
els (1.5-2 nM) that are used in our standard assay conditions. 
Since the estimated IC50 value was in the same concentration 
range as the enzyme, the inhibition constant may prove to be even 
lower. When the enzyme concentration was lowered or raised, the 
IC50 value changed as well (Figure ID). This suggests that we 
can give an estimate for the apparent K, value to be in the range 
of 100 pM to 1 nM. The exact value of such a constant is diffi- 
cult to determine using available methods. This extremely tight 
inhibition complex should exhibit a very low off-rate, which is 
an important factor in modeling the Sicl phosphorylation and 
degradation switch, as discussed below. 

Sicl IS EFFICIENTLY PH0SPH0RYLATED WITHIN THE Sic1/Clb5-Cdk1 
COMPLEX 

The dynamics of the switch are also shaped by the rates of 
Sicl phosphorylation by Clnl,2-Cdkl and Clb5,6-Cdkl. We 
have previously shown that both Clb5-Cdkl and Cln2-Cdkl 
efficiently phosphorylate the N-terminal part of the truncated 



non-inhibitory version of Sicl (Sicl AC) (Koivomagi et al., 
201 la). Substrate specificity was shown to be gained by an inter- 
action between a hydrophobic docking site in Clb5 and RXL 
motifs in Sicl, and in the case of Cln2-Cdkl via a novel hydropho- 
bic LLPP-docking motif in Sicl (Koivomagi et al., 201 lb; Bhaduri 
and Pryciak, 2011; Koivomagi and Loog, 2011). Furthermore, by 
detailed mapping of the paths of the multisite phosphorylation 
cascades we found that the Clb5 specificity toward different sites, 
including the pair S76/S80, which form the critical diphospho- 
degrons at the end of the cascade, was much more pronounced 
compared to Cln2. These observations led to the conclusion that 
although Cln2-Cdkl can phosphorylate Sicl, it cannot trigger 
Sicl degradation without the help of Clb5-Cdkl. This conclu- 
sion was supported by in vivo experiments showing that, in the 
absence of B-type cyclin (Clb) activity, the Clns alone were not 
able to cause the degradation of Sicl. Thus, Clb-Cdkl activity is 
the driving force of Sicl degradation at the Gl/S transition. This 
underscores the likely importance of a double-negative feedback 
mechanism whereby Clb5-Cdkl phosphorylates Sicl, leading to 
its destruction and thereby releasing further active Clb5-Cdkl 
from the inhibitory Sicl-Clb5-Cdkl complex. 

This model requires that Clb5-Cdkl can efficiently phospho- 
rylate Sicl bound to the inhibited Clb5-Cdkl complex. To test 
this, we performed a Sicl phosphorylation assay using partly 
inhibited Clb5-Cdkl — that is, we tested the ability of Clb5 
to phosphorylate a Sicl-Clb5-Cdkl complex (Figures 2 A-D). 
Additionally, the phosphorylation of non-inhibitory Sicl AC was 
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FIGURE 2 | Sicl is efficiently phosphorylated within the 
Sic1/Clb5-Cdk1 complex, when an excess of Clb5-Cdk1 is added. 

Assays with Cln2-Cdk1 complexes are included for comparison. (A) 
Autoradiography showing the phosphorylation of purified Sidwt and 
SidAC by Clb5-Cdk1 and Cln2-Cdk1 complexes using 15 nM Sid and 



30 nM cyclin-Cdkl complexes. (B) The activity of Clb5-Cdk1 was 
determined using 2.5 ilM histone H1, with and without 15nM Sidwt in 
the assay mixture. (C,D) Quantification of the experiment presented in 
panel (A). For comparison, phosphorylation rates were normalized using 
cyclin-Cdkl activity units obtained in panel (B). 
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assayed in the same conditions. In these assays, 30 nM Clb5- 
Cdkl or Cln2-Cdkl was mixed with 15 nM Sicl. To quantify the 
exact fraction of uninhibited Clb5-Cdkl activity, a control sub- 
strate histone HI was included in a separate assay (Figure 2B). 
Clb5-Cdkl phosphorylated both forms of Sicl efficiently and at 
a considerably higher rate than Cln2-Cdkl. The rate of phospho- 
rylation of the full length Sicl was similar to that of Sicl AC for 
both Cln2 and Clb5. These data indicate that the N-terminus of 
Sicl is freely accessible for phosphorylation when Sicl is bound 
in the inhibitory complex, and that emerging free Clb5-Cdkl, 
released from the inhibitory complex as Sicl is degraded, could 
accelerate Sicl phosphorylation and degradation, thus creating a 
double-negative feedback. 

We have also observed that a low level of intramolecular 
Sicl phosphorylation takes place within the inhibitory Sicl-Clb5- 
Cdkl complex during its formation (data not shown). It is par- 
tially also the likely reason why the full length Sicl showed slightly 
faster phosphorylation rate compared to Sicl AC (Figures 2C,D). 
Indeed, structural modeling studies do not exclude the possibil- 
ity that the active site of Cdkl may be partly accessible during 
the inhibitor binding (Barberis et al., 2005a,b; Barberis, 2012a,b). 
The possible role of these phosphorylation events in the switch is 
currently under investigation. 

THE MODEL OF THE CDK/CKI-CONTROLLED G1/S SWITCH 

To quantitatively analyze the impact of the inhibition strength 
and the phosphorylation rates on the abruptness and the degree 
of irreversibility of the Gl/S switch, we constructed a mini- 
mal model of Sicl phosphorylation and degradation (Figure 3, 
Table 1). As a first step, we analyzed the steady state output of 
the model. This is important to understand the system in gen- 
eral and to analyze the transition thresholds for the conditions 



(e.g., slow growth, the arrest, or stress) when the degree of irre- 
versibility of the switch may become especially important. As a 
second step we analyze the behavior of the switch in the normal 
dynamic circumstances of the cell cycle. 

To simplify our model and avoid equation systems that lack 
steady state solutions, we did not include the multisite phos- 
phorylation steps, but considered the net phosphorylation and 
degradation process as a single step. This was based on our 
experimental data showing that accumulation of multiply phos- 
phorylated species is a rapid, semi-processive process in vitro, and 
that multiply phosphorylated forms of Sicl do not accumulate 
in vivo after cells are released from an a-factor arrest (Koivomagi 
et al, 2011a). The latter observation suggests that Sicl degrada- 
tion is rapid compared to phosphorylation. Therefore, we can 
consider the rate of phosphorylation and degradation to be a 
single value that is mostly determined by phosphorylation rates. 
Also, due to the negligible accumulation of multi-phosphorylated 
forms of Sicl and prompt degradation in response to increasing 
Cdkl activity in vivo, we omitted the counteracting phosphatase 
from the model. This does not mean that the phosphatase does 
not take part in the process but that, for simplicity, its effect 
was included in the net rate constant of the whole process. The 
first of the three ordinary differential equations of the model 
describe free Sicl levels, including the association and dissocia- 
tion of the inhibitory complex and Cln2-dependent phosphory- 
lation/degradation (Table 1). The second equation describes the 
dynamics of free Clb5 concentration, which is controlled by the 



Table 1 | Equations describing the minimal model of the G1/S switch. 
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y 2 = [Clb5 - Cdkl total] - Y3 

/i = [Sid total] - Y3 

Kinetic constants (the values presented below are those used in 

Figure 4. Parameter values that are varied in other figures are given in 

the legends or on the figure panels): 

fci = 0.01 min -1 ; dissociation rate of the inhibitory complex 

l<2 = 0.1 nM -1 x min -1 ; association rate of the inhibitory complex 

k 3 = 0.001 nM -1 x min -1 ; the net rate constant for Clb5-dependent 

phosphorylation and degradation of Sid 

k A = 0.00001 nM -1 x min -1 ; the net rate constant for Cln2-dependent 
phosphorylation and degradation of Sid 

= 10nM x min -1 ; the basal synthesis rate of Sid 
fc 6 = 0.01 min -1 ; the basal degradation rate of Sid 

The description of the system is presented in Figure 3. For the phase dia- 
grams presented in Figures 5-7, the steady-state solutions for the depen- 
dencies between Sicl total and Clb5-Cdkl total and Clb5-Clb5-Cdkl free and 
Clb5-Cdk1 total were derived. 
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FIGURE 3 | Minimal model of the Gl/S switch. Parameters used for 
diagrams and simulations presented in the figures below are listed in 
Tables 1 and 2. Black crosses designate the phospho-dependent 
degradation of Sid via the SCF-proteasome pathway. Gray crosses 
designate basal degradation of Sicl. 
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association-dissociation rates of the inhibitory complex and also 
by the fast release of Clb5-Cdkl from the inhibitory complex 
upon Sicl destruction (Verma et al., 2001). The third equation 
describes the concentration changes of the inhibitory complex. 
The conservation equations define the total levels of Clb5 and 
Sicl. Finally, the model includes slow basal synthesis and degrada- 
tion of Sicl . This was based on the finding that although the major 
Swi5-induced transcription of SIC1 peaks at the exit from mitosis 
(Schwob et al., 1994), SIC1 transcription is still observed through- 
out the entire cell cycle (Knapp et al., 1996; Aerne et al., 1998). 
Taken together, although considerably simplified, this minimal 
model contains the basic and universal core elements of a switch 
containing a double-negative feedback between a CKI and a CDK. 
Additionally, the system contains a second CDK complex (the 
equivalent of Cln2-Cdkl) that is not inhibited by CKI but acts 
as an initiator of the switch. 

The values for the phosphorylation rate constants were derived 
from our experimentally determined rates of the phosphoryla- 
tion of Sicl and other specific model substrates by cyclin-Cdkl 
complexes (Loog and Morgan, 2005; Koivomagi et al, 2011a,b). 
A Km value of 3.3 u,M and a k cat value of 73 min -1 were deter- 
mined for Clb5-Cdkl -dependent phosphorylation for Sicl (data 
not shown). Considering that there are about six to seven phos- 
phorylation steps, and only a slight phosphatase counteraction, 
the net rate constant for the process of Clb5-dependent phospho- 
rylation (kcibs) and degradation of Sicl was taken to be 0.001 
nM _1 x min -1 . Although Cln2-Cdkl showed only a 3-4-fold 
lower rate for net phosphorylation of Sicl in the Sicl-Clb5-Cdkl 
complex (Figure 2), it was shown to have a considerably weaker 
ability to phosphorylate the priming site T33, and also, to finalize 
the phosphorylation cascade, thus creating a di-phosphodegron 
pS76/pS80 (Koivomagi et al., 2011a). Therefore, the net rate con- 
stant for Cln2-dependent phosphorylation (kci n 2) was initially 
taken to be 0.00001 nM _1 x min -1 , two orders of magnitude 
lower than the corresponding value for Clb5-Cdkl. The impor- 
tance of this difference on the properties of the switch will be 



further analyzed below. Note, as explained above, that these rate 
constant are the net constants for both phosphorylation and 
degradation. 

The parameters defining the inhibitory interaction between 
Sicl and Clb5-Cdkl were derived from the inhibition experiment 
presented in Figure 1. Assuming that the diffusion-limited associ- 
ation rate of an average protein-protein interaction that involves 
an intrinsically disordered binding partner is in the range of 10 7 
to 10 8 M _1 x min , then a Kd of a few hundred pM would 
yield a dissociation rate for the Sicl-Clb5-Cdkl inhibitory com- 
plex in the range of 10~ 3 to 10~ 2 min -1 . The effect of varying 
this parameter will be studied in the context of the switch. The 
basal synthesis and degradation rates were taken to be relatively 
slow, and for the phase diagrams below these values were chosen 
to give a steady-state concentration of Sicl of 1 u,M in the absence 
of cyclin-Cdkl complexes. 

When the steady-state solution of this ODE system is presented 
as a plot of total Clb5 levels vs. the fraction of free Clb5 at a fixed 
Cln2 concentration, a steep hysteresis loop can be observed, indi- 
cating that the system is a bistable switch (Figure 4A). The levels 
of total Sicl vs. total Clb5 also exhibit hysteresis (Figure 4B). 
These diagrams show the two states of the switch separated by two 
saddle node bifurcations. The state with high total Sicl and low 
free Clb5-Cdkl corresponds to Gl, and the state with low total 
Sicl and high levels of free Clb5-Cdkl corresponds to S-phase. If 
the cell is in Gl and the total Clb5 concentration is gradually ris- 
ing, the phase paths, as shown by the arrows, reach a threshold 
upon which there is a jump in free Clb5-Cdkl steady-state value 
and a drop in Sicl levels. This is the on-switch for S-phase. This 
bistable system sets a high threshold for Clb5, thereby filtering 
noise, before a sharp transition that ensures an abrupt and coher- 
ent start of S-phase processes. Also, the switch between the two 
states is highly irreversible and resistant to transient fluctuations 
in Clb5 levels. Once activated, it is hard to reverse the switch, as it 
would take a several-fold drop in Clb5 concentration before Sicl 
can re-accumulate to switch the system back to the Gl state. 




FIGURE 4 | Phase diagrams showing the testability of the system at 
steady state. Parameter values are provided in Table 1. The dark gray arrows 
show possible paths for a gradual increase of free Clb5-Cdk1 while entering 
the S-phase, and possible reversible paths (light gray arrows) in the event 
of stochastic decreases in Clb5 levels. (A) Phase diagram showing the 
dependence of free Clb5-Cdk1 activity on total Clb5-Cdk1 levels. When a 



threshold is reached, the system switches itself into the state corresponding 
to S-phase, as indicated by the upward pointing arrow. (B) Phase diagram for 
the same system with the same parameter values as in panel (A), showing the 
dependence of total Sid concentration on total Clb5-Cdk1. The threshold level 
for a drop of Sid corresponds to the same threshold level of total Clb5-Cdk1 
that triggers the jump in steady state levels of free Clb5-Cdk1 in pane! (A). 
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TIGHT INHIBITION OF Clb5-Cdk1 BY Sid AND EFFICIENT 
PHOSPHORYLATION OF Sic1-Clb5-Cdl<1 BY FREE Clb5-Cdk1 ARE 
CRUCIAL FOR SETTING A DISCREET THRESHOLD FOR THE G1/S 
TRANSITION 

To understand the basic features of the switch and to determine 
the optimal balance of crucial parameters that are required to 
maintain the bistable nature of the system, we studied the effect of 
inhibition strength and phosphorylation rates on the behavior of 
the phase diagrams. In the diagrams presented in Figures 5A-D, 
inhibition strength was varied by introducing different values for 
the dissociation rate of the inhibitory complex. As the dissocia- 
tion rate increased, the bistable behavior of the system gradually 
weakened and the switch evolved to a smooth transition that can 
easily slide back and forth (Figure 5D). This kind of a behavior 
would potentially be very harmful for a transition like Gl/S: only 
a slight backward deviation of the total Clb5 levels would reverse 
the cell to Gl, interrupting processes like the initiation of repli- 
cation, and removing the re-replication block. Thus, our model 
predicts that between a K, of 10-100 nM the on-threshold con- 
tinues to decrease, and eventually the system loses its hysteresis, as 
well as its switch-like response (Figures 5C,D). In this range, Sicl 
becomes a readily overcome barrier for the accumulating Clb5. 

Next, in the diagrams presented in Figure 6, we explored 
the effect of relative phosphorylation rates. To understand the 
impact of the experimentally determined strikingly weaker ability 
of Cln2-Cdkl to phosphorylate the critical output degron sites 
in Sicl, compared to Clb5-Cdkl (Koivomagi et al., 2011a), we 
first varied the Cln2-dependent phosphorylation rate relative to a 
fixed Clb5-dependent rate (Figures 6A-C). Importantly, steeper 
hysteresis was observed when the Cln2-dependent phosphoryla- 
tion rate was considerably lower compared to the rate of Clb5. 
Increasing the relative Cln2-dependent rate resulted in less hys- 
teresis, lowered initial Sicl levels, and lowered the Clb5 threshold 
for the transition (Figures 6A,B). Thus, the timing of the transi- 
tion could potentially be tuned by Cln2 activity. However, when 
Cln2 specificity toward Sicl approaches the specificity of Clb5, the 
requirement for Clb5 signal almost ceases to exist (Figure 6C). 
Similarly, as in the case of a weak K, (above), the threshold of 
Total (Clb5) decreased and the discontinuity of the Clb5-Cdkl 
switch was lost. Thus, both tight inhibition and cyclin specificity 
are important factors for an abrupt transition and a high degree of 
irreversibility. Interestingly, the effect of Cln2 specificity on bista- 
bility was more pronounced compared to the effect of an increase 
in K, of the same order. 

Thus, in the system described in Figure 3, where both Cln2 
and Clb5 are able to send Sicl for degradation, the signal strength 
of Cln2 can potentially modulate timing by lowering the thresh- 
old of Clb5 required to trigger the switch. However, our recent 
results indicate that Cln2-Cdkl cannot phosphorylate the last 
phospho-degron sites in Sicl at a sufficient rate to cause degra- 
dation in vivo. On the other hand, Cln2-Cdkl can efficiently 
phosphorylate the initial residues in the Sicl cascade, and this 
primes Sicl for efficient subsequent phosphorylation by Clb5- 
Cdkl, which will finalize the cascade to send Sicl to degradation 
(Koivomagi et al., 2011a). To study how such a priming effect 
could influence the switch, we gradually raised the phosphoryla- 
tion rate constant for Clb5 in the absence of Cln2 activity. In this 
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FIGURE 5 | (A-D) The steady state phase diagrams showing the effect of 
the inhibition strength on the bistability of the Gl/S switch. The dissociation 
rates of the inhibitory complex were increased in 10-fold increments from 
0.01 to 10min _l as indicated. In the model, these values correspond to K, 
values of 0.1-100 nM. Arrows indicate the paths of the system at the G1/S 
transition. 
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FIGURE 6 | (A-C) Phase diagrams showing the effect of changes in the 
relative ability of Cln2-Cdk1 to directly cause the degradation of Sid. The 
rate constant values for Cln2-Cdk1 were gradually increased as indicated. 



simplified way, the different constants would roughly mimic dif- 
ferent levels of Cln2 as a priming kinase (higher rate constants 
being equivalent to primed Sicl). The diagrams in Figure 7 reveal 
that, similarly to the scenario in Figure 6, priming by Cln2 would 
influence the threshold of the switch. However, the abruptness of 
the transition is not lost within a realistic range of rate constant 
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FIGURE 7 | Phase diagrams showing the potential effect of 
Cln2-dependent priming phosphorylation on the bistability of the 
system. The rate constant value of 0.001 nMxmin -1 for Clb5-Cdk1 
exemplifies a system with no Cln2-dependent priming while increasing 
values of this constant mimic increased levels of Cln2-dependent priming 
effect. (A) Phase diagrams showing the dependence of free Clb5-Cdk1 
activity on total Clb5-Cdk1 levels. (B) Phase diagrams showing the 
dependence of total Sid concentration on total Clb5-Cdk1. The color 
coding is the same for both panels. 



values for Clb5-Cdkl. Such tuning of the threshold of the switch 
by priming can act as a mutual control system between the 
Gl-CDK and S-CDK activities. If there is not enough S-phase 
CDK activity, the Sicl degradation switch is not triggered. If there 
are sufficient levels of the S-phase CDK activity, but for some rea- 
son the levels of Gl-CDK drop, the transition would be delayed. 
Thus, a priming mechanism has the advantage of connecting the 
Cln and Clb signals, allowing tuning of decision-making by both 
pathways. When Cln and Clb signals are independent, there is 
a possibility that the switch will be triggered before sufficient 
amounts of Clb5 are accumulated. In the priming model, Sicl 
degradation and lifting the barrier from S-phase always requires 
sufficient levels of Clb5 or other early expressed B-type cyclins. 
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CYCLIN CONCENTRATIONS AND SYNTHESIS RATES 

Although the phase diagrams provide a good general charac- 
terization of an ODE model describing a bistable system, they 
don't necessarily fully reflect the way the system behaves in reality. 
Some rates could be too slow relative to the others, in which case 
steady states would not be established. In the present case, the 
synthesis of Clb5 is fast, with a peak accumulation time of only 
around 10-20 min (Koivomagi et al., 2011b). Therefore, there 
is a possibility that relatively slow phosphorylation rates or dis- 
sociation rates of the inhibitory complex would not allow the 
formation of the steady states described by the phase diagrams. To 
address this concern, we investigated the dynamics of our model. 
We numerically simulated a modified form of the ODE system 
(Table 2). We searched the literature to supply realistic synthe- 
sis rates and peak concentrations of the cyclins and Sicl. We 
obtained experimentally estimated endogenous concentrations of 
Clnl, Cln2, Cln3, Clb5, Clb6, and Sicl. The most detailed studies 
have been performed by western blotting of ProtA-tagged pro- 
teins and by using quantitative calibration of the blotting signals 
with ProtA standards (Cross et al, 2002). Second, a study pre- 
senting a global analysis of protein abundances in yeast provides 
quantifications based on the western blotting of TAP- tagged pro- 
teins (Ghaemmaghami et al., 2003). Using the ProtA-tag method, 
the estimated average number of Sicl molecules per diploid cell 



Table 2 | Equation system used for numerical simulation of the time 
courses at G1/S in Figures 8A-C, 9A-D, and 10A. 



— = -k 2 yi y 2 - fc^i Ya + fci Y3 
dt 

dy2 

—r = k e - k 2 y\ n + k 3 y 2 y 3 + k 4 y 3 y 4 + ki y 3 
-jT = k 2 yi(y 6 - y 3 ) - k 3 y 2 y 3 - k A y 3 y A - k[ y 3 



— = -ktyi y A - k 3 y 2 y 3 - k 4 y 3 y 4 

dy 6 . 

-dF = ks 

Yi = [Sid FREE ] 

K2 = [Clb5-Cdk1 F REE] 

y 3 = [Sic1/Clb5-Cdk1] 
y 4 = [Cln2-Cdkl] 
y 5 = [Sid TO tal] 
y 6 = [Clb5- Cdk1 total] 

The kinetic constants (values presented below are those used in 
Figure 8A. Parameter values for other figures are given in the legends or 
on the figure panels): 

kf = 0.01 min -1 ; dissociation rate of the inhibitory complex 
k 2 = 0.1 nM -1 x min -1 ; association rate of the inhibitory complex 
k 3 = 0.001 nM -1 x min -1 ; the net rate constant for Clb5-dependent 
phosphorylation and degradation of Sid 

fct = 0.00001 nM -1 x min -1 ; the net rate constant for Cln2-dependent 
phosphorylation and degradation of Sid 
kj = 100 nM x min -1 ; Cln2 synthesis rate 
^8 = 30 nM x min -1 ; Clb5 synthesis rate 

The initial values for y1 (SIcIfree) and y5 (Si'cItotal) in two simulations 
were 1 fiM. 



in asynchronous cultures was reported to be 214. The study using 
the TAP-tag in haploid cells provided an estimate of 768 copies of 
Sicl per cell. To calculate the realistic initial concentration value 
for Sicl for the ODE simulations, one needs to take into account 
the predominantly nuclear localization of Sicl (Huh et al., 2003). 
The nuclear volume of haploid and diploid cells is thought to 
be 2 and 4 fl, respectively (Cross et al, 2002). Thus, the average 
estimated nuclear concentration of Sicl in asynchronous cultures 
would be 363 nM. Since Sicl is present during only a fraction 
of the cell cycle, we considered that a reasonable peak concen- 
tration would be three times that of the average abundance in 
asynchronous cells. Therefore, a realistic concentration of Sicl at 
its peak in Gl would be around 1 |xM. In fact, this value falls into 
the same range as the steady state value of Sicl in the absence of 
Clb5 and Cln2 in the phase diagrams presented in Figure 7. 

Analogous calculations for Clb5, whose localization is also 
predominantly nuclear (Shirayama et al., 1999), give peak con- 
centration of 1138nM (115nM for Clb6). The cyclins Clnl and 
Cln2 are localized both in cytoplasm and in nucleus (Miller and 
Cross, 2000; Edgington and Futcher, 2001; Landry et al, 2012), 
and therefore the cell volume for calculating their concentration 
was taken to be 100 fl for the ProtA study (diploid cells) and 50 fl 
for the TAP study (haploid cells). Interestingly, the obtained peak 
values for Clnl and Cln2 were about 164 and 33 nM, respectively, 
which is much lower than the peak values for Clb5 and Clb6. This 
fact would additionally support the model that Clnl, 2 may have a 
less critical role in catalyzing Sicl degradation compared to rela- 
tively abundant Clb5. Cln3 levels were recorded only by the study 
using ProtA, and their concentration in Gl, taking into account 
the nuclear localization, was in the range of 269 nM. 

Based on these calculations and the half-times of accumula- 
tion of Cln2 and Clb5, the concentrations of Cln2 and Clb5 in this 
model were set to rise at 30 and 100 nM per minute, respectively. 
To obtain a clear and formal picture of the effect of accumulation 
rates on the dynamics of the switch, a constant linear accumu- 
lation was introduced (Table 2). During the short time window 
of cyclin accumulation, the basal synthesis and degradation rates 
of Sicl used in the ODE system for the phase diagrams would be 
too slow to significantly affect the dynamics of Sicl degradation. 
Therefore, these rates were excluded to lower the complexity of 
the model. 

TIGHT INHIBITION AND EFFICIENT Clb5-DEPENDENT Sid 
PHOSPHORYLATION IS IMPORTANT FOR THE ABRUPTNESS 
OF THE SWITCH UNDER REALISTIC DYNAMIC CIRCUMSTANCES 

The simulation of the ODE system presented in Figure 8A repre- 
sents a situation with tight inhibition and high kcib5/kcin2 ratio, 
analogous to the one presented in Figure 4. Although exhibit- 
ing a time course with a sigmoidal shape, a fairly slow transition 
from low to high Clb5-Cdkl activity level is observed. Thus, 
although the phase diagrams can provide a general description 
of the system, the transitions between the alternative states may 
happen less abruptly. We also investigated the effect of lowering 
the kcib5/kcin2 ratio. Similarly to the conclusion made based on 
the phase diagrams above, a relatively faster Cln2-dependent rate 
decreased the ultrasensitivity of the signal response and slowed 
the transition (Figure 8B). Also, analogously to the tendency 
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FIGURE 8 | Time course simulations to analyze the abruptness of the 
transition. The equation system and parameter values are provided in 
Table 2. (A) The simulation with the basic set of key parameters for 
inhibition and phosphorylation, similar to those used in the phase diagrams 
in Figure 4. (B) The effect of faster Cln2-dependent phosphorylation rate 



on the switch. The rate constant value for Cln2-Cdk1 was taken to be 
equal to the "k" value of Clb5-dependent phosphorylation. (C) The effect 
of weaker inhibition strength on the switch. The dissociation rate of the 
inhibitory complex was raised to 10min _l , which in the model 
corresponds to a K,- value of 100 nM. 



observed in Figure 5, increasing the K, value in these dynamic 
circumstances reduced the abruptness of the switch (Figure 8C). 

If Cln2-Cdkl was left out of the model and Clb5-Cdkl had 
its usual rate constant value of 0.001 nM -1 x min -1 , the sys- 
tem behaved like the one in Figure 8A (kci n 2 = 0.0001 nM -1 x 
min -1 ), with free Clb5-Cdkl only starting to appear when the 
concentration of total Clb5 surpassed the initial concentration 
of Sicl at 1000 nM (Figure 9A). In both of these cases the start 
of Sicl degradation corresponds to the point at which Clb5 



surpasses stoichiometry with Sicl . The slow Sicl degradation pro- 
file, in this case, shows that the kcibs value of 0.001 nM -1 x 
min -1 alone would not be sufficient to cause the transition to 
be abrupt. However, in such a situation S-phase would still be ini- 
tiated with only a few minutes delay. This is in agreement with 
our experimental observations that the mutation of the Clnl,2- 
specific docking motif in Sicl (the LLPP motif) introduced about 
a 5-10 min delay in Sicl degradation after the release of the cells 
from the Gl arrest (Koivomagi et al., 2011a). 
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FIGURE 9 | Time course simulations showing the potential effect of 
Cln2-dependent priming phosphorylation on the abruptness of the 
switch. The equation system and parameter values are provided in Table 2. 
(A) A system with no priming effect in the absence of Cln2. The basal rate 



constant value of 0.001 nlVT 1 x min -1 for Clb5-Cdk1 was used. (B-D) 
Different priming effects were mimicked by using the following "k" values for 



Clb5-Cdk1 
min" 1 (D). 



0.01 nlvr 1 x min" 1 (B), 0.03 nM" 1 x min" 1 (C). and 0.1 nrVT 



Frontiers in Physiology | Systems Biology 



December 2012 | Volume 3 | Article 459 | 10 



Venta et al. 



The G1/S switch 



When the effect of Cln2 was introduced in the capacity of 
a primer for Clb5, it was possible to tune the threshold of 
the transition (Figures 9A-C). In these simulations, different 
phosphorylation constant values for Clb5 present examples of 
different priming effects created by Cln2. Importantly, Cln2- 
dependent priming introduced abruptness to the switch. It is 
plausible that priming (phosphorylation of sites by Cln2-Cdkl 
that Clb5-Cdkl can subsequently utilize as docking sites for Cksl- 
dependent finalization of the cascade) creates a mixture of Sicl 
forms whose net phosphorylation rate by Clb5 is optimal to 
cause a rapid decline of Sicl levels. Indeed, lowering the num- 
ber of steps in the cascade that Clb5 must pass through would 
effectively raise the net rate of the Clb5-dependent process. Also, 
priming of the non-Cksl -dependent slow initial step by another 
kinase would lessen the burden on Clb5. Interestingly, when the 
kcibs value was raised further to O.lnM -1 x min -1 , abrupt- 
ness was retained, but the Clb5 threshold for the Gl/S transition 
would be about three times lower than the Clb5 peak value 
(Figure 9D). This scenario would create a weak noise barrier 
and would also lack the potential for a burst activation of S- 
phase with the full potential of accumulating Clb5. This situation 



would potentially cause genomic instability, a phenotype char- 
acteristic of the SIC1A strain (Lengronne and Schwob, 2002). 
Thus, the estimated range of net Clb5-dependent phosphoryla- 
tion rates of 0.01 nM min -1 for the mixture of primed forms 
of Sicl is in an ideal window to achieve abruptness while main- 
taining a noise-filtering threshold. It is conceivable that multisite 
phosphorylation systems provide wide combinatorial flexibility 
for tuning and optimization of these kinetic properties during 
evolution. 

We also tested the situation in which Cln2-Cdkl alone was the 
kinase responsible for degradation of Sicl, and Clb5-dependent 
positive feedback was omitted (kcibs = 0). The abruptness of 
the switch was considerably reduced and the whole transition 
followed a path of a slow and gradual change (Figure 10A). 
This strongly suggests that the optimal mechanistic basis and the 
source of the ultrasensitive signal response of the switch would be 
the double-negative feedback effect of Clb5. If Clnl,2 were to act 
alone, the transition would not be sufficiently abrupt and there is 
no guarantee that it would not happen prematurely with insuffi- 
cient Clb5 signal to promote S-phase. However, an earlier model 
of Sicl phosphorylation and degradation based the responsibility 
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FIGURE 10 | Time course simulations for systems where only Cln2 is 
responsible for phosphorylation and degradation of Sicl. (A) A model 
similar to that in Figures 8 and 9 was used, except that the rate constant 
value for Clb5-Cdk1 toward Sid was taken to be zero. The rate constant for 



Cln2-Cdk1 was 0.001 nM -1 x min -1 , which is the same value used for 
Clb5-Cdk1 in Figure 8A. (B) A system that requires Cln2-dependent multisite 
phosphorylation for degradation of Sicl. The sequential phosphorylation of six 
sites in Sid was taken to be the output signal for degradation of Sicl. 



www.frontiersin.org 



December 2012 | Volume 3 | Article 459 | 11 



Venta et al. 



The G1/S switch 



of Sicl phosphorylation entirely on Cln2 (Nash et al, 2001; Verma 
et al., 1997a). The ultrasensitive feature of the switch in this 
model was later built on a theoretical presumption of equal rates 
of distributive phosphorylation and the delayed accumulation 
of multiply phosphorylated forms of Sicl (Deshaies and Ferrell, 
2001). To test if this scheme would improve the abruptness of 
the switch, we introduced six equally Cln2-specific phosphory- 
lation steps into the model, with phosphorylation at all six being 
required for degradation. Clb5-dependent positive feedback effect 
was excluded as in Figure 10A (for other details of this model 
see Table 3). The simulation of this system did indeed create a 
longer lag for the onset of the transition, but it did not improve 



the abruptness of the response (Figure 10B). In reality, the path 
of the transition would look even less switch-like because the N- 
terminus of Sicl contains several crucial phosphorylation sites 
bearing suboptimal CDK consensus motifs and exhibits low indi- 
vidual site specificity with Cln2-Cdkl (Koivomagi et al., 2011a). 
These simulations predict that although the delayed accumula- 
tion of multiply phosphorylated forms could create a threshold 
for Sicl degradation, multisite phosphorylation would not add 
considerably to the abruptness of the switch. This conclusion is 
in agreement with the results of a theoretical analysis predict- 
ing that multisite phosphorylation may provide a threshold for 
a transition but is generally a poor switch (Gunawardena, 2005). 



[Clb5-Cdk1] 



Table 3 | Equation system used for simulation of the time course for Cln2-dependent multisite phosphorylation of Sid presented in Figure 10B. 
Scheme 

[Cln2-Cdk1] [Cln2-Cdk1] [Cln2-Cdk1] [Cln2-Cdk1] [Cln2-Cdk1] /C/„, 

[Sid] ► [1P-Sid] ► [1P-Sic1] ► [1P-Sid] ► [1P-Sid] ► [5P-Sid] ~^-c 9 . 

k1 k2 k3 k4 k5 k6"» 

kon I kon I kon I kon I kon I kon I r ^K„S., 

*H H -H ~H -M -H ^ [6 K 1] 

|Cln2-Cdk1] [Cln2-Cdk1] [Cln2-Cdk1] [Cln2-Cdk1] [Cln2-Cdk1l *$r**<< 

[Sid ► [1P-Sid ► [2P-Sid ► [3P-Sid ► [4P-Sid ► [5P-Sid 

Clb5-Cdk1] k1 Clb5-Cdk1] k2 Clb5-Cdk1] k3 Clb5-Cdk1] k4 Clb5-Cdk1] k5 Clb5-Cdk1] [Clb5-Cdk1] 

The kinetic constants 

kon = 0.1 ; — the association rate of the inhibitory complex. To lower the complexity, the dissociation rate was considered very slow and not included in 
this version of the model. 

kl = 0.001 ; k 2 = 0.001 ; k 3 = 0.001 ; k 4 = 0.001 ; and k 5 = 0.001 ;— the rate constants for Cln2-dependent phosphorylation steps. 

k§ — 0.001 ; — the rate constant for Cln2-dependent phosphorylation of 5P-Sic1 and 5P-Sic1/Clb5-Cdk1 including the net rate for the subsequent fast 

SCF-mediated degradation. 

/cNi = 30; fcSi = 100; the synthesis rates of Cln2 and Clb5, respectively. 
The variables 

K(i)i K(2)i Y{3)', Y(4)', Y(5)', K(6>; — Sid forms bearing 0, 1, 2 or 5 phosphates, respectively. The form with 6 phosphates was degraded by SCF and due 

to the fast rate of degradation this form was not expressed as a separate variable. To lower the complexity, only a sequential 6-step process was 
considered in the model. 

VO)', V&y, V&)', Voo)', YmY' K(12); — different Sic1/Clb5-Cdk1 inhibitory complexes involving Sid with 0, 1, 2 or 5 phosphates, respectively. 

y<i3)) y<i4)) YOS)', — concentrations of free Clb5-Cdk1 , Cln2-Cdk1, and total Clb5-Cdk1, respectively. 

The initial value for 0P-Sic1 (yp)) was taken to be 1, the initial values for the rest of the variables were taken to be zero. 

The ODE system 

dy m /dt = -ki x y (U ) x y m - kon x (y (15) - y (7) - y (8) - y (9) - y m - y (m - y (12) ) x y m 
dy^/dt = fc| x y (14) x y m - k 2 x y m x y {2) - kon x (y (15) - y a) - y (8) - y (9) - y (10) - - y (12) ) x Y(l) 
dy (3) /dt = k 2 x y (14) x y (2) - k 3 x y (14) x y {3) - kon x (y (15) - y (7) - y (8) - y (9) - y (10) - K<ii) - /(12)) x V(3) 
dy m /dt = k 3 x y (14) x y (3) - k 4 x y (14) x y w - kon x (y- (15) - y a) - y (S) - y (9) - K(10) - K<ii) - /(12)) x V(4) 
dy^/dt = fc 4 x y (14) x |/(4) - k s x y (14) x y {5) - kon x (y- (15) - y a) - y (S) - y (g) - y (W) - y m - y 02) ) x y {5) 
dyw/dt = k s x y (14) x y ( 5) - k 6 x y (14) x y {6) - kon x (y- (15) - y a) - y (8) - y {3) - y m - y m - y (12) ) x / (6) 

dy (7) /dt = kon x (y (15) - y a) - y [S) - y (S) - y 00) - y m - y m ) x y m - k\ x y (yl) x (y aS) - /(i3) - /(8) - V(9) - V(10) - Yav) - V(12)) 

dy^/dt - kon x (/(is) - y (7} - y m - y (9) - y m - y m - y m ) x y (2) + (k : x y 04) ) x (y (15) - y (13) - y m - y (9) - y m - y m - y a2) ) - (k 2 x y (14) ) x 

(K(15) - K(13) - Y(7) - K(9) - V(10) - Ym) - V(12)) 

dy<9)/dt = kon x CVciS) - Ym ~ Ym ~ YO) ~ Yoo) - Y(n) - V(12)) x K(3) + (fe x y 04) ) x (y (15) - y (13) - y m - y (9) - y m - y m - y a2) ) - (k 3 x y (14) ) x 
(K(15) - K(13) - K(8) - YP) - Y(W) - K(H) - K(12)) 

dy { iO)/dt = kon x (v (15) - y (7) - y (8) - y (3) - y m - y m - y m ) x y m + (k 3 x v (14) x (y (15) - y 03) - y (B) - y {7) - y m - y m - v, 12) ) - (k 4 x y (14) ) x 
(V(15) - K(13) - K(8) - K(9) - Yd) - Veil) - K(12)) 

dy m /dt = kon x (y 0S) - y (7) - y (3) - y (a) - y m - y m - y (m ) x y (S) + (k 4 x y (U )) x (K(15) - Y()3) - Y(S) - Y&) - Ym - - K(12)) - (h x K(14)) x 

(K(15) - V(13) - K(8) - K(9) - K(10) - Y(7) - K(12)) 

dya 2) /dt = kon x (y (16) - y 0) - y m - y (9) - y 00) - y m - y m ) x y {6) + (k 5 x y (14) ) x (y 05) - y (13) - y (B) - y (9} - y m - y g) - y m ) - (fe x y ( i 4) +) x 
(K(15) - K(13) - K<8) - K(9) - K(io> - Ym) - VO)) 

dy 03) /dt = kB^ - kon x (y (15) - y (7) - y {3 ) - y (9 ) - Y(W) - Veil) - K<12)) x (y m + y {2) + y (3) + y (4) + y (5) + y {6) ) + k 6 x y (14) x y 02) 
dy (U )/dt= kNi 
dyc\5)/dt= kB } 
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DISCUSSION 

The double-negative feedback model of the Gl/S switch in bud- 
ding yeast involves two CDK complexes and a CDK inhibitor 
Sicl that, at the same time, is the substrate of CDK. We aimed 
to explore which properties such a mutual antagonism would 
require to generate a useful switch. Based on mathematical mod- 
eling and phosphorylation assays we found that, for the transition 
to exhibit the best switch-like qualities, it is important to have 
tight inhibition of Clb5-Cdkl by Sicl and efficient phosphory- 
lation of the Sicl-Clb5-Cdkl inhibitory complex. The low K; 
ensures that a tight and nearly stoichiometric barrier sets a 
threshold for the Clb5 signal. The efficient phosphorylation and 
degradation rates ensure that the double-negative feedback would 
gain the momentum required to create a switch-like increase in 
free Clb5-Cdkl in response to total Clb5. We also studied the rel- 
ative impact and different roles of Cln2 and Clb5 in the switch. 
In this respect we analyzed several possible scenarios: first, a sce- 
nario in which both Cln2 and Clb5 are capable of independently 
triggering the degradation of Sicl; second, a case in which Cln2 
phosphorylates Sicl to prime further Cksl -dependent phospho- 
rylation by Clb5 but is not able to cause the degradation of 
Sicl by itself; third, we modeled the situations in which either 
Cln2 or Clb5 alone was responsible for phosphorylation and 
degradation of Sicl; finally, we explored the impact of multisite 
phosphorylation on the switch. The general conclusion derived 
from these case studies was that the optimal role for Cln2-Cdkl 
in the switch would not be to act as a kinase activity responsible 
for abrupt degradation of CKI, but rather to act as a prim- 
ing signal that triggers the feedback loop for the emerging free 
Clb5-Cdkl. 

Clb5-dependent positive feedback would be predicted to play 
a role in Sicl degradation at the Gl/S transition, since there is 
expected to be a point when the accumulating fraction of free 
Clb5-Cdkl will overcome and phosphorylate Sicl molecules that 
are bound to the inhibitory complex. This feedback is clearly 
expected to have some contribution to the abruptness of the 
switch. However, until recently, the extent of this effect had not 
been quantitatively estimated. In addition to our recent find- 
ings, Tang and colleagues used live cell fluorescent microscopy 
to show that Clb5,6-dependent positive feedback promotes the 
switch-like properties of the Gl/S transition (Tang, pers. commu- 
nication). Additionally, recent work by Barberis and colleagues 
has shown the involvement of Clb5 and also later B-type cyclins 
in degradation of Sicl (Barberis et al, 2012). The authors of this 
work proposed a model where a surviving fraction of Sicl could 
be involved even post-Gl/S in the regulation of transcriptional 
waves of later cyclin genes. In general, these findings further argue 
against a solely Cln2 -dependent mechanism of Sicl degradation 
and emphasize the relatively major importance of B-type cyclins 
in this process. 

The mutual control between Clb5 and Cln2 in tuning the 
threshold point and preventing a premature Gl/S transition 
emphasizes the advantage of Cln2 as a primer rather than as a 
qualitatively similar kinase signal relative to Clb5. In addition, 
we have found that there could be other kinases besides Clnl,2 
involved in priming of Sicl for Clb5, since some phosphorylation 
of an N-terminal site is seen in Gl-arrested cells (Koivomagi et al., 



201 la). The inability of Clnl,2 to degrade Sicl in vivo would likely 
also be partly caused by the much lower nuclear concentration of 
Cln2 compared to Clb5. Thus, it is possible that Clnl,2-Cdkls 
play only a minor role in the switch. 

The initial idea that Clnl,2-Cdkl alone could be responsible 
for degradation of Sicl was mainly based on an experiment in 
which about 2-fold overexpression of the non-degradable form 
of Sicl, relative to the endogenous levels of Sicl, did not inhibit 
endogenous Sicl degradation but was able to block the start of 
replication (Verma et al., 1997a). However, our recent experiment 
showed that about 4-fold overexpression (our unpublished quan- 
tifications) of the non-degradable inhibitory domain of Sicl fully 
stabilized endogenous Sicl (Koivomagi et al., 2011a). Thus, it 
seems that at relatively low levels of the inhibitor the active Clb5- 
Cdkl is able to accumulate to the levels that exceed the threshold 
capable of degrading Sicl, but not to the levels required to ini- 
tiate S-phase. Another study reported that when non-degradable 
Sicl accumulates to about the same levels as endogenous wild- 
type Sicl, cells are still capable of initiating DNA replication with 
only a minor delay (Cross et al, 2007). It seems that if the bar- 
rier set by the non-degradable inhibitor is at the same levels as 
wild-type Sicl, then Clb5 will accumulate and can still trigger 
S-phase. If the barrier is about 2-fold, it can degrade Sicl but 
cannot start S-phase. This raises the interesting possibility that 
the Clb5-Cdkl threshold level for initiation of Sicl degradation 
is lower than the threshold for triggering S-phase by phospho- 
rylation of Sld2 and Sld3. This would ensure that the transition 
always starts before the initiation of replication. Interestingly, 
a similar, highly resolved substrate phosphorylation order has 
been demonstrated for the accumulation of mitotic cyclin-CDK 
activity in both budding yeast (Oikonomou and Cross, 2011) 
and human cells (Deibler and Kirschner, 2010; Gavet and Pines, 
2010). 

The minimal model introduced in this study could serve 
as a core system that could help to understand Gl/S control 
modules in other organisms. The crucial CKI at the Gl/S tran- 
sition in human cells, the protein p27Kipl (Hengst et al., 1994; 
Polyak et al, 1994; Toyoshima and Hunter, 1994), must be 
phosphorylated by CDK to be recognized by the SCF ubiqui- 
tin ligase for subsequent degradation by the proteasome pathway 
(Pagano et al, 1995). It is believed that a positive feedback 
effect is created by CyclinE-Cdk2-dependent phosphorylation 
of the phosphodegron residue at Thrl87 (Sheaff et al., 1997). 
Furthermore, the phosphorylation of Thrl87 is facilitated by 
priming phosphorylation at Tyr88 by tyrosine kinases of the 
extracellular mitogen initiated signaling pathways (Chu et al, 
2007; Grimmler et al., 2007; Jakel et al., 2012). The priming 
grants CDK the access to Thrl87 and partly activates the inhib- 
ited CyclinE/Cdk2 by loosening the inhibitory interaction. In 
fact, human tyrosine kinases, by transmitting growth factor or 
cytokine signals, and Clnl,2-Cdkl in yeast, by transmitting the 
signal of the cell-size control system, both play similar roles in the 
cell cycle by initiating S-phase. In both of these cases, the priming 
activity is not inhibited by CKI. Thus, it is possible that the varia- 
tions of the proposed priming mechanism coupled to the mutual 
negative feedback of the CKI/CDK system could be function- 
ally conserved across species. Therefore, the modeling approach 
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presented here to estimate kinetic parameters that influence the 
abruptness and irreversibility of the Gl/S switch may provide 
general insight into Gl/S control in humans. This knowledge is 
especially important since perturbations in the regulatory mod- 
ule governing the Gl/S is widely observed in cancer (Chu et al, 
2008). 

METHODS 

PROTEIN PURIFICATION 

The TAP method was applied for purification of cyclin-Cdkl 
complexes as described previously for Clb5-TAP-Cdkl (Puig 
et al, 2001; Ubersax et al, 2003). 3HA-Cln2-Cdkl was puri- 
fied according to published protocols (McCusker et al, 2007) 
using rabbit polyclonal antibody against the haemagglutinin epi- 
tope (Labas). N-terminal 6His-tagged recombinant Siclwt and 
SiclAC proteins were purified by cobalt affinity chromatogra- 
phy. Cksl was purified as described in Reynard et al. (2000). 
Protein concentrations were measured by colloidal coomassie 
G-250 using BSA as a standard. 

KINASE ASSAYS 

The general composition of the kinase assay mixture was as fol- 
lows: 50 mM HEPES pH 7.4, 5 mM MgCl 2 , 150 mM NaCl, 0.1% 
NP-40, 20 mM imidazole, 2% glycerol, 2 mM EGTA, 0.2 mg/ml" 1 
BSA, 500 nM Cksl, and 500 |xM ATP [with added [y- 32 P]ATP 
(Perkin Elmer)]. For the inhibition assay twelve reactions with 
Siclwt concentrations ranging from 10 to 0.04 nM and two reac- 
tions containing no inhibitor were used. As a reference substrate, 
bovine histone HI (Upstate) with final concentration 2.5 (jlM 
was used. The concentration of kinase complex in the assay was 
1.5-2 nM. Before starting the reaction enzyme-inhibitor com- 
plex was formed (2min). Reaction was started by adding HI. 
Two timepoints, 7 and 14 min, were collected and analyzed. For 
the phosphorylation assays presented in Figure 2, 15 nM Siclwt 



or SiclAC and 30 nM kinase complex were used. Four time- 
points were collected: 20s, 40s, 60s, and 80s. For demonstrating 
enzyme activity and Sicl inhibition ratios 2.5 uM HI was used as 
a substrate. For data analysis GraphPad Prism 5.0 software was 
used. 

MODELING 

The models were constructed based on mass action kinetics 
using values of the kinetic constants derived from biochemi- 
cal experiments. Because Sicl concentrations in the model were 
several times below the Km values for Cdkl -dependent phos- 
phorylation of Sicl, velocities were expressed using kcat/R-M "is a 
second order rate constant (vo = k cat /KM[S][E]) instead of the 
Michaelis-Menten equation. The kinetic schemes of the models 
are presented in Figure 3 (for model simulations in Figures 4-7) 
and Table 3 (for model simulations in Figure 10B). The time 
course simulations in Figures 8, 9, and 10A were based on the 
scheme in Figure 3 with the exceptions that the synthesis and 
the basal degradation rates of Sicl were omitted and linear syn- 
thesis rates for Cln2 and Clb5 were added. For the steady state 
diagrams we solved the equation system using MATLAB and the 
steady state values for plotting the diagrams were calculated with 
a step of 0.01 of the total Clb5 levels. The time course simula- 
tions for Figures 4-10 were performed using the MATLAB ode45 
solver. 
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